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The electron-phonon matrix elements relevant to the Raman intensity and Kohn anomaly of 
the G band are calculated by taking into account the effect of the edge of graphene. The analysis 
of the pseudospin reveals that the longitudinal optical phonon mode undergoes a strong Kohn 
anomaly for both the armchair and zigzag edges, and that only the longitudinal (transverse) 
optical phonon mode is a Raman active mode near the armchair (zigzag) edge. The Raman 
intensity is enhanced when the polarization of the incident laser light is parallel (perpendicular) 
to the armchair (zigzag) edge. This asymmetry between the armchair and zigzag edges is useful 
in identifying the orientation of the edge of graphene. 
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interaction, pseudospin, gauge field 



O 

(N 
> 
CO 
OS 

in 



o 



X 



1. Introduction 

Graphene is a unique material since its electron mo- 
tion is governed by a special equation similar to the rel- 
ativistic massless Dirac equation, while a nonrelativistic 
equation is common in condensed matter physics. 1, 2 -* The 
electron motion is modified by the electron-phonon (el- 
ph) and electron-light interactions, which are fundamen- 
tal issues in discussing the transport, 1,2 ) electronic, 3 ' and 
optical properties 4 ' 5 -* of graphene. The goal of this paper 
is to show that an asymmetry of the Raman spectra for 
r point longitudinal and transverse optical phonon (LO 
and TO) modes, both of which are known as the Raman 
G band, appears near the edge of graphene. There are 
two fundamental orientations for the edge of graphene, 
zigzag and armchair edges, and a general edge shape is 
considered to be a mixture of them. 6, 7 > The asymme- 
try is useful in identifying the orientation of the edge of 
graphene by Raman spectroscopy. 

In Raman spectroscopy, we irradiate laser light onto 
a sample and observe the intensity of the inelastically 
scattered light. The energy difference between the in- 
cident laser and the inelastically scattered light corre- 
sponds to the energy of a Raman active phonon mode 
due to the energy conservation. The el-ph interaction is 
essential for the Raman process. Further, the el-ph inter- 
action can modify the energy and life-time of the phonon 
mode, which is known as the Kohn anomaly. 8 ' Evidence 
for Kohn anomalies is found in the phonon dispersion of 
carbon nanotube, 9 - 1 graphene, 10,11 ' and graphite. 12 - 1 By 
examining the Kohn anomaly for the G band of carbon 
nanotube, 13 ' a feature of the el-ph interaction such as 
the chirality dependence of the el-ph interaction upon 
the Kohn anomaly has been clarified. 14 - 1 In this paper, 
we calculate the el-ph matrix elements relevant to the 
Raman intensity and Kohn anomaly of the G band of 
graphene within effective-mass approximation by includ- 
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ing the effects of the edge of graphene and polarization 
direction of an incident laser (and a scattered) light. 

This paper is organized as follows. In § 2, we show the 
Hamiltonian including the el-ph interaction with respect 
to the r point optical phonon modes and the electron- 
light interaction. In § 3 and § 4, we calculate the matrix 
elements for the el-ph and electron-light interactions by 
taking into account of the presence of the zigzag and 
armchair edges, respectively. The self-energy of the LO 
mode is estimated in § 5 and the phonon self-energy for 
general edge shape is discussed. Finally, we propose two 
models representing the electronic states at the interior 
of a graphene sample and calculate the self-energies for 
those models in § 6. In § 7, we discuss the relationship 
between our result and experimental results, and sum- 
marize the results. 

2. Hamiltonian 

Let v I / K(r) [^k' (r)] be the wave function for an electron 
near the K [K'] point, the energy eigen equation for an 
electron near the Fermi energy of graphene is written as 



H 



* K (r) 
*K'(r) 



E 



* K (r) 
*K'(r) 



(1) 



The wave function ^>k(^) [^K'( r )] is two-component 
structure, which results from that the hexagonal unit 
cell contains two carbon atoms [A atom (•) and B atom 
(o) in Fig. 1]. The total Hamiltonian H including the el- 
ph interaction with respect to the T point LO and TO 
modes, and the electron-light interaction is given by 15 - 1 



H = v F 



(p 



- eA) 







(p-Ao 



(2) 



Here vp is the Fermi velocity, momentum operator p = 
-ih(d x ,d y ), a = (a x ,a y ) and a' = (—a x ,a y ) where a x , 
<j y and a z are Pauli matrices. We take x and y axes 
as shown by the inset in Fig. 1(a). The electromagnetic 
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gauge field A enters into the Hamiltonian through the 
substitution p — > p eA where — e is the charge of 
electron. A uniform field A can represent the incident 
laser light and the scattered light in the Raman process. 
The el-ph interaction is represented by the deformation- 
induced gauge field A q = (A 1 ^, A^). 15 ^ It can be shown 
that A^ and A^ are expressed in terms of a change of the 
nearest-neighbor hopping integral from an average value 
- 7o , <5 7 o, a , as 16 - 18 ) 



V F A% = (57o,l - - (#70,2 + 570,3) : 

v f a v = "^r (^0,2 - ^70,3) • 



(3) 



Here a (= 1,2,3) for S-f , a denotes the direction of the 
bond (sec the inset of Fig. 1), and <57o, a is caused by 
atomic displacements by the T point optical phonon 
modes. Note that A q is uniform for the F point q = 
phonons, while A q depends on the position r as A q (r) 
for phonons with q 7^ 0. 14 ' Although an additional 
deformation-induced gauge field due to a local modu- 
lation of the hopping integral originating from a defect 
appears in a realistic situation, we ignore it in eq. (2) for 
simplicity. 

3. Zigzag Edge 

First, we calculate the matrix element relevant to the 
Raman intensity near the zigzag edge. The scattering or 
reflection of an electron at the zigzag edge is intravalley 
scattering, 19 ' and therefore we can consider the K and 
K' points separately. Let us consider the electrons near 
the K point. The Hamiltonian is given by 



H K = v F a • (p + A q - eA) . 



(4) 



We specify the deformation-induced gauge field A q for 
the LO and TO modes near the zigzag edge. The vi- 
brations of carbon atoms corresponding to the T point 
LO and TO modes are shown in Figs. 1(a) and 1(b). By 
assuming that the perturbation £70, a is proportional to 
the change in the bond length, we have #70,1 = and 
^70. 2 = ~<57o. 3 for the LO mode, while #70,1 = — 2<57o.2 
and 570, 2 = ^70,3 for the TO mode. Using eq. (3), we see 
that A q for the LO mode is written as A^ Q = (0,^4 q ) 
with v F A^ = a/3^70,2, while A q for the TO mode is 
written as A^ Q = (A q ,0) with v F A^ = — 3<57o,2- Note 
that the direction of A q for the LO (TO) mode is per- 
pendicular (parallel) to the zigzag edge. The direction 
of A q is perpendicular to the direction of atom dis- 
placement. 9 ' 2 °) Thus, the el-ph interaction in eq. (4), 
#lo/to = V F (T ■ Al / X0 , is rewritten as 



(5) 



H% = v F A«a v , 
H z T l = v F A«a x , 

for the LO and TO modes, respectively. 

The el-ph matrix element is given as the expectation 
value of the el-ph interaction with respect to the energy 
eigenstate for the unperturbed Hamiltonian, — v F cr ■ 
p. The energy eigenstate with wave vector k in the con- 
duction energy band is written in terms of the plane wave 



(a) LO mode 



A LO 



(b) TO mode 




^70,2 = -#70,3 

x <57o,i = 



£70,2 = 670,3 
670,1 = -2670,2 



Fig. 1. The displacement vectors for the LO and TO modes are 
shown in (a) and (b), respectively. The displacement vectors of 
the LO (TO) mode are parallel (perpendicular) to the zigzag 
edge. The direction of the deformation-induced gauge field A q 
is perpendicular to the direction of atom displacement. 



e ik r and the Bloch function $ k as $ k (r) = Ae ikr $ k , 
where A is a normalization constant satisfying N 2 V = 1, 
V is the area (volume) of the system, and 



(6) 



Here 0(k) is the angle between the vector k and the k x - 
axis (see Fig. 2). The expectation values of a x , a y , and 
a z with respect to define the pseudospin. Since a x = 
<* c k KI*k} = cos#(k), a v = = sin0(k), and 

°z = (^kl a z|^k) — 0' the direction of the pseudospin of 



-k' 



(a x ,a y ,a z ) = (cos0(k),sin0(k),O), 



(7) 



is within the (k x ,k y ) plane and parallel to the vector k 
(see Fig. 2). Owing to the presence of the zigzag edge par- 
allel to the x-axis, the wave function near the zigzag edge 
is a standing wave given by a sum of the incident wave 
$ k (r) and the reflected wave $ k /(r) with k' = (k x , —k y ) 
as 



*£(r) = ^(*£(r) + *£,(r)) 



(8) 



Strictly speaking, it is necessary to add the relative phase 
between $ k (r) and $ k /(r) in order that ^/ k (r) may sat- 
isfy the boundary condition for the zigzag edge. However, 
this phase gives no contribution to the matrix elements 
of interest in the present investigation, and therefore we 
omit it. Note that the normalization of eq. (8) is adopted 
for k y 7^ 0. Some complications arise when k v = 0. For 
example, when k y — and k x < 0, localized wave func- 
tions of edge states 21 ' should be used, which is explained 
in Appendix A. 

The el-ph matrix element from a state in the conduc- 
tion band to the same state is given as (* k |# L o/ TO | , I' k )- 
Using eqs. (5) and (8), the pseudospin, and 0(k') = 
— 0(k), we obtain 

(*kl^T§l*k)=«FA q cosc?(k). (10) 

This result shows that the Raman intensity of the LO 
mode is negligible compared with that of the TO mode 
at zigzag edges. For eq. (9), (^kl^s/l^k) can be rewritten 
as a sum of two components, ($ k |<7 y |$ k ) -I- ( ( f > k 'l cr i,| < I' k ') 1 
since cross terms such as ($ k |<T y |$ k/ ) vanish. Because 
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Fig. 2. The zigzag edge reflects the wave vector k = (k x ,k y ) to 
k' = (k x , —k y ), and two wave functions 'f^(r) and ( r ) form 
a standing wave. Note that the y-component of the pscudospin 
is flipped by the zigzag edge. 



Here, we have used |\I>£) = er z | , I'£), a y a z — ia x , a x cr z — 
—i<Ty, and eq. (11). It is noteworthy that it is mainly 
the electrons near the fc x -axis [9 (Is.) « or tt] that can 
participate in the Raman process taking place near the 
zigzag edge since both the el-ph matrix element [cq. (10)] 
and the optical transition amplitude [cq. (12)] are pro- 
portional to cos#(k). Let us define the angle between 
the laser polarization and the zigzag edge as 9 (see the 
inset in Fig. 3), then A y = |A|sin0 and the Raman 
intensity is proportional to |M opt (A)| 2 cx sin 2 9. The 
9-dependence of the square of the optical transition am- 
plitude is plotted as the dashed curve in Fig. 3. 



sin0(k) and ($ k ,K|$ k ,) = sin0(k') = 



— sin0(k), the y-component of the pseudospin for the 
incident wave ^(r) is reflected as shown in Fig. 2. Thus, 
we have 



o. 



(11) 



due to the cancellation of the y-component of the pseu- 
dospin between the incident k-state and the reflected k'- 
state. Similarly, in eq. (10), (^kl^l^k)' i s written as 
a sum of two components, ($£|o- x |$£) + (<&£., |<r x | <&£.,). 
Because (3> k |<7 x |$ k ) = cos#(k) and (3> k /|CT x |$ k /) = 
cos6>(k') = cos#(k), we obtain eq. (10). 

The Raman intensity depends on the polarization 
of the incident laser light 22 ) and on that of the scat- 
tered light. The electron- light interaction is given by 
iJg- m = —vpecr ■ A in cq. (4). The optical absorption 
occurs with amplitude M p*(A) = (* k |-fflTl*k>> where 
^(r) is the wave function in the valence energy band, 
which is related to * k (r) via *k( r ) = CT «*k( r )- 0n 
the other hand, the optical emission occurs with ampli- 
tude (^kl^K m l*k)' which is simply the complex con- 
jugate of M opt (A). Thus, the polarization dependences 
of the incident and scattered light are the same. Here, 
let us examine the polarization dependence of the inci- 
dent light. The direction of A = (A x ,A y ) corresponds 
to the direction of the polarization of the electric field. 
The polarization of the incident laser light should be 
perpendicular to the zigzag edge within a graphene 
plane, i.e., A^ = (0,^4^), in order to populate photoex- 
cited electrons effectively. This argument follows from 
M opt (Aj_) = -iv F e A y cos 0(k) for A± = (0,A y ), while 
M p*(A||) = for A,, = (A x ,0) because 

M^(A ± )^-v F eA y (^l\a y m) 

= -v F eA y (^\a y a z \^) 

= -iv F eA v (^\a x \^) 

= -iv F eAy cos6>(k), (12) 



and 



M°P t (A||) = - WF eA x (^|a x |^) 

= -v F eA x {^l\a x a z \^l) 

= iv F eA x (^l\a y \^l) 

= 0. (13) 



3 
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zigzag only 

armchair only 

(pure : r = 0) 

zigzag@armchair 

(partial : r = 0.5) 

zigzag@armchair 

(random : r = 1) 



Laser Polarization 

.A 




6(degree) 



Fig. 3. The polarization dependence of the square of the optical 
transition amplitude (|M opt (A)| 2 ) is plotted as a function of the 
angle of laser polarization (©) with respect to the orientation 
of the edge. For a pure zigzag (armchair) edge, the intensity is 
maximum when the laser polarization is perpendicular (parallel) 
to the edge. "zigzag@armchair" denotes the case when zigzag 
edges are introduced into part of a perfect armchair edge. We 
have used eq. (29) with r = (armchair only), r = 0.5 (partial), 
and r = 1 (random: a mixture of zigzag and armchair edges). 



The Kohn anomaly is relevant to the el-ph ma- 
trix element for electron- hole pair creation, i.e., 
<*kl#u)/Tol*k>- Usin S *k(r) = *** k (r), we rewrite 
the matrix element as (*kl^Lo/TO CTz l*k)- F rom c 1- (5), 
we have 



H i!o a z = iv F A^a x , 



H TO a z = - iv FA x <J y , 



(14) 



where a x a z 



-la 



y and OyO z 



have thus shown that 



■ ia x have been used. We 
is proportional to a y as 
well as that H™q is proportional to a y . From eq. (11), we 
see that the TO mode is unable to transfer an electron 
in the valence band into the conduction band, that is, 
the TO mode does not decay into an electron-hole pair, 
and therefore the Kohn anomaly for the TO mode is 
negligible compared with that for the LO mode. 

4. Armchair Edge 

Next, we calculate the matrix element relevant to the 
Raman intensity near the armchair edge. Suppose that 
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the armchair edge is located along the y-axis, then the 
armchair edge reflects an electron with k = [k x , k v ) near 
the K point into the state with k' = (k' x , k' y ) = (—k x , k y ) 
near the K' point, where k and k' are measured from the 
K and K' points, respectively. The negative sign in front 
of k x for k' is due to the momentum conservation. One 
may consider that an intervalley process is unconnected 
with the r point LO and TO phonons. Note, however, 
that we should consider the K and K' points simultane- 
ously in the case of an armchair edge since the reflection 
of an electron by the armchair edge is an intervalley scat- 
tering process as shown in Fig. 4. 



k'„ 



(a) LO mode 



(b) TO mode 





^70,2 = #70,3 
670,1 = -2^70,2 



670,2 = -670,3 
57o,i = 





Pseudospin 






jjn 


e'(k') ] , 


\ K 


k x { K.' 


T — J 




Fig. 5. The displacement vectors for the LO and TO modes are 
shown in (a) and (b), respectively. The displacement vectors of 
the LO (TO) mode are parallel (perpendicular) to the armchair 
edge. The direction of A q is perpendicular to the direction of 
atom displacement. 




state near the K' point can be expressed as 



1 
71 



1 

_ e -i0'(k' 



(18) 



Fig. 4. The armchair edge reflects the wave vector k = (k x ,k y ) 
of one valley into k' = (—k X; k y ) of another valley, and the two 
wave functions of the different valleys form a standing wave. 
The pseudospin is unchanged by the armchair edge. Note that 
the pseudospin for states near the K' point is not parallel to the 
vector k', while the pseudospin for states near the K point is 
parallel to the vector k. 



We specify the deformation-induced gauge field A q for 
the LO and TO modes near the armchair edge. The vi- 
brations of carbon atoms for the T point LO and TO 
modes are shown in Fig. 5. We have 670,1 = —2670,2 
and 670,2 = 670,3 for the LO mode, while 670,1 = and 
670,2 = —670,3 for the TO mode. Using eq. (3), we see 
that A q for the LO mode is written as A^ Q = (A q ,0) 
with vpA^ — —3670.2, while A q for the TO mode is writ- 
ten as A^, = (0, A q ) with v F A^ = ^670,2- Thus, from 
eq. (2), we see that the el-ph interaction 



armchairedge where 8'(k') is defined through k' x + ik' — \k'\e l9 ( k ). 



Since the armchair edge reflects the state with k — 
(k x ,k y ) into the state with k' = (—k x ,ky), we have 
the relation 6>'(k') = vr - 0(k) (see Fig. 4). By substi- 
tuting this into eq. (18), we see that the Bloch func- 
tion of eq. (18) becomes of eq. (6), which explains 
cq. (17). The pseudospin for the eigenstate near the 
K' point is given by ($£., Icr^ | ) = — cos6>'(k') and 
($£.' I <7 y |<f>£./) = sin6»'(k'). Thus, the pseudospin for the 
K' point is not parallel to the vector k', as shown in 
Fig. 4, although the pseudospin for states near the K 
point is parallel to the vector k. Using #'(k') = it — 9(k), 
one can see that the pseudospin is preserved under the 
reflection at the armchair edge (see Fig. 4). 

Using eqs. (16) and (17), it is straightforward to check 
that 



<^|ffr O m |*k>=«F^COS0(k), 



(19) 



rrarm 

-"lo/to 



vf 



■ A q 

-^■LO/TO 






A q 

LO/TO 



(15) 



is rewritten as 



rrarm 



v?A* 



= v F A« 










— (7 



This result shows that the Raman intensity of the TO 
mode is negligible compared with that of the LO mode. 
The absence of the Raman intensity of the TO mode 
results from the interference between two valleys, namely, 
the opposite signs in front of a y for the K and K' points 
of Hfg 1 in eq. (16). 

The interaction between the light and the electronic 
states is given by 



(16) 



H cm (A) = -w F e 



a ■ A 







r'-A 



(20) 



for the LO and TO modes, respectively. 

The wave function is given by a sum of the plane wave 
at the K point and the reflected wave at the K' point as 

*£(r) = -=,(*^_ ifc .J. (17) 



from eq. (2). The optical absorption amplitude is given 
by M°p*(A) 



(*^|F m (A)|l v k ), where 



V2 



*k(r) 







*£(r). 



(21) 



Note that the Bloch function is the same ($£) f° r both 
the K and K' points. In fact, the Bloch function for a 



If the polarization of the incident laser light is per- 
pendicular to the armchair edge A^ = (A x ,0), then 
XAjJI'I^) vanishes owing to the cancellation 
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between the K and K' points. The polarization of 
the incident laser should be parallel to the armchair 
edge, i.e., A|| = (0,A y ), in order to populate photoex- 
cited electrons effectively because (*£.|-ff em (A||)|\I> k ) = 
— iv F eA y cos#(k). Note that it is mainly the electrons 
near the fc^-axis [0(k) w or n] that can participate the 
Raman process taking place near the armchair edge since 
both the el-ph matrix element [eq. (19)] and the opti- 
cal transition amplitude are proportional to cos#(k). By 
defining the angle between the laser polarization and the 
armchair edge by O (see the inset in Fig. 3), we have 
A y = |A|cos0, and we see that the Raman intensity 
is proportional to |M opt (A)| 2 oc cos 2 0. The polariza- 
tion dependence of the Raman intensity for the armchair 
edge is opposite that for the zigzag edge, as shown in 
Fig. 3, from which the orientation of the edge may be 
determined experimentally. 

The el-ph matrix element for the Kohn anomaly is 
given by (^kl^Lo/TO^I^k)- From e 1- ( 16 )' we nave 







H^a z = iv F A y 



(22) 



It has thus been shown that the TO mode does not un- 
dergo a Kohn anomaly because the matrix element van- 
ishes owing to the sign difference between the K and K' 
points with respect to a x . 

5. Energy Difference Between LO and TO 
Modes 

In this section we calculate the energy difference be- 
tween the LO and TO modes. The renormalized phonon 
energy is written as a sum of the unrenormalized energy 
haj and the self-energy. Since the TO mode does not un- 
dergo a Kohn anomaly, the self-energy of the TO mode 
vanishes. Thus, the energy difference between the LO and 
TO modes is the self-energy of the LO mode, which is 
given by time-dependent second-order perturbation the- 
ory as 



l(* k |g L ar o m l*k)l 2 
hu + E? h + iS 



X (/h - /e) • 



(23) 



where the factor of 2 originates from the spin degeneracy, 
/e,h = (l+cxp((£; ' h -£; F )/fcBT)~ 1 is the Fermi distribu- 
tion function, E F is the Fermi energy, 8 is a positive in- 
finitesimal, E° (E h ) is the energy of an electron (a hole), 
and Eg 1 = E^-E^ = 2hv F \k\ (> 0) is the energy of an 
electron-hole pair. Note that the summation index ^ k 
in cq. (23) is not restricted to only interband (E eh ^ 0) 
processes but also includes intraband (E ch — 0) pro- 
cesses. Thus, the self-energy can be decomposed into two 
parts, U(u,E F ) = n intor (w, E F ) + n intra (w, E F ), where 
n intcr (u;, E F ) includes only interband electron- hole pair 
creation processes satisfying E ch ^ 0. 

In the adiabatic limit, i.e., when u = and <5 = 
in eq. (23), by substituting eq. (19) into cq. (23), it is 



straightforward to show that, at T = 0, 

V / A q \ 2 

n intra (o,£ F ) = -^ f-^J \e f \, 

V / A q \ 2 
n intCT (o,i? F ) = -^ f-^J (E C -\E F \) 



(24) 



where E c is a cutoff energy. Note that IT ntra (0, E F ) does 
not vanish because (/h — f e )/E^ 1 ^ in the limit of 
££ h -> 0, while in the nonadiabatic case, n mtra (cj, E F ) 
vanishes since (/h — f c )/fou = in this limit. It is only 
the interband process that contributes to the self-energy 
in the nonadiabatic case. Lazzeri and Mauri 11 -* pointed 
out that 11(0, E F ) does not depend on E F in the adiabatic 
limit owing to the cancellation between n intra (0, E F ) and 
n intcr (0 ,E F ). This shows that the adiabatic approxima- 
tion is not appropriate for discussing the E F dependence 
of the self-energy. In the nonadiabatic case, at T = 0, it 
is a straightforward calculation to obtain (see Appendix 
B for derivation) 

Re[n(w,£ F )] = 



V (A% 
it \ h 



E r 



|£ F |-^ln 



\E F \ - !f 



\E, 



(25) 



The Fermi energy dependence is given by the last two 
terms. 10, n ) The first term is linear with respect to E F 
and the second term produces a singularity at E F = 
±tko/2. These terms express the nonadiabatic effects. 23 ^ 
Recently, Saitta et al. 24 ^ have pointed out that large 
nonadiabatic effects are found to be more ubiquitous in 
layered metals such as CaCg and MgB 2 . 
For the case of E F = 0, eq. (25) becomes 

V /A q \ 2 

Re[n(w,0)] = --M-jfJ E c . (26) 

The self-energy depends on the cutoff energy E c . The 
value of E c cannot be determined within the effective- 
mass model. We assume that E c is of the order of half of 
the 7r bandwidth (10 eV); see §7 for a detailed discussion 
of the value of E c . Using the harmonic approximation 
for the displacement of the carbon atoms, 25 ) we obtain 
y/TQA%/h\ w 2 x lCr 2 !" 1 (see Appendix B), where N u 
denotes the number of hexagonal unit cells. Since V can 
be written as N U S where S is the area of a hexagonal 
unit cell, we obtain Re [II (w, 0)] w —6 mcV. Thus, the 
difference in the Raman shift between the (Raman ac- 
tive) TO mode near the zigzag edge and the (Raman ac- 
tive) LO mode near the armchair edge is approximately 
50 cm -1 . In a realistic system, the actual magnitude of 
the self-energy may be much smaller than this value. For 
example, a typical edge is a mixture of zigzag and arm- 
chair edges, for which the energy difference between the 
LO and TO modes is lower. 

Here, let us introduce zigzag edges into part of a per- 
fect armchair edge at x = and examine the effect of 
the randomness of the edge shape on the Raman inten- 
sity and phonon self-energies. Then the standing wave 
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near the rough edge is approximated by 



+ ae 



ik'- 



+ ze 



ik'-: 



(27) 



where k' = (—k x ,k y ) and |a| 2 + |z| 2 = 1. The wave 
function ^(r) reproduces cq. (17) for the case when 
(a, z) = (1,0). Note that |z| 2 /|a| 2 (= r < 1) can bc con- 
sidered phenomenologically as the ratio of the number of 
zigzag edges to that of armchair edges in the rough edge, 
and r = 1 [(\a\, \z\) — (l/\/2, l/\/2)] represents the case 
that armchair and zigzag edges are equally distributed 
along the y-axis. It is a straightforward calculation to 
obtain 



1 



(*kl^b m |*k) = ^cos0(k) T 



(n\H^\n) = v F A*sm9(k)- 



'■ + \zf 

+ M 2 



(28) 



+ a 2 + z 



2 ' 



These matrix elements show that the self-energy for 
the LO mode becomes Re [II (w, 0)] /4 for the case of 
r = 1. On the other hand, the self-energy of the TO 
mode, which is zero for the case of r = 0, becomes 
Re [II(u;,0)] /4 for the case of r = 1. The differences in 
the Kohn anomalies for the LO and TO modes disappear 
for the case of r = 1. Moreover, the Raman intensity of 
the TO mode increases, while the Raman intensity of 
the LO mode decreases. As a result, the G band exhibits 
a single peak. The intensity of the G band is given as 
the sum of the LO and TO modes. Since the intensity of 
each mode is four times smaller than that of the LO mode 
near the pure armchair edge, the total intensity of the G 
band should be two times smaller than the Raman inten- 
sity near the pure armchair edge. Note that for a general 
value of (a, z), the energy difference between the LO and 
TO modes is given by (|a| 2 — |z| 2 )Rc [II(w,0)]. It is also 
a straightforward calculation to obtain the polarization 
dependence of the optical transition amplitude, 



|M opt (A)| : 



cos 2 O r 2 sin 2 O 



oc 



(29) 



(1 + r) 2 (1 + r) 2 ' 

This dependence is plotted for two cases, r = 1 and r = 
0.5, in Fig. 3. 

6. Bulk and Edge 

In the case of an infinite periodic graphene system 
without an edge, the self-energies of the LO and TO 
modes are the same and given by Re [II (w, Ep)] in 
eq. (25). Moreover, no asymmetry between the LO and 
TO modes in the Raman intensity is expected. The rea- 
son why the LO and TO modes do not exhibit any differ- 
ence in Raman spectra is that graphene is a homo-polar 
crystal with two atoms per unit cell, and hence there is 
no polar mode, similar to the case of Si. Thus, the LO 
and TO modes are degenerate and contribute equally to 
the single peak of the G band (see "Periodic" in Fig. 6). 
Note that a slight change in the spring force constant due 
to a uniaxial strain applied to a graphene sample can re- 
solve the degeneracy between the LO and TO modes. In 
this case, the unrenormalized energy hw for the LO mode 



is not identical to that for the TO mode. However, even 
for this case, we can expect that the self-energies and 
Raman intensities for the LO and TO modes are similar 
to each other. Thus, we can see two peaks for the LO 
and TO modes with similar intensity, as was observed 
by Mohiuddin et al. 2 ^ 

Since an actual sample is always surrounded by an 
edge, it is interesting to consider whether or not the in- 
terior of a graphene sample can be considered as an in- 
finite periodic graphene system without the edge. If the 
wave function in the interior region is given by a su- 
perposition of the incident and reflected states, then it 
is reasonable to assume that the wave function is ap- 
proximated by eq. (27) with (|a|,|z|) = {l/y/2, l/v/2), 
since it is probable that the edge is a random mixture 
of zigzag and armchair edges. The peak positions of the 
LO and TO modes in the Raman shift are indicated by 
"Random" in Fig. 6. We speculate that the peak position 
for an actual sample appears between the peaks labeled 
"Periodic" and "Random" . An estimation of the effective 
distance from the edge at which the effect of interference 
on the pseudospin discussed so far can survive will be a 
subject of further investigation. 



— O- 



Zigzagedge ; L0 



Bulk 



"It 



Random 

-V- 



Periodic 



Armchair edge 

hw + Re [n(<j, ())] 



Raman Shift 



hw + -Re [no, 0)] 



Fig. 6. The horizontal lines indicate the Raman shift for the case 
of Ep = 0. (top and bottom lines) The Raman peak taken near 
the zigzag (armchair) edge appears only for the TO (LO) mode 
indicated by the solid circle. The peak for the TO mode does not 
accompany the broadening because the TO mode decouples from 
the electron- hole pairs, (middle line) The Raman peak appears 
at hui -f Re [TI(uj, 0)] in the case of an infinite periodic graphene 
system without an edge ("Periodic"). If the effect of the electron 
reflection at the edge survives in the interior of a graphene sample 
("Random"), the Raman peak is expected to appear at tiui + 
Rern(«,0)] /4. 



7. Discussion and Conclusions 

Here, we discuss the relationship between our result 
and experimental results. Cangado et al. observed that 
the Raman intensity of the G band for a nanoribbon 
has a strong dependence on the incident light polariza- 
tion. 27 ) They showed that the Raman intensity is max- 
imum when the polarization is parallel to the edge of a 
nanoribbon. Their result is consistent with our result for 
the armchair edge, but not consistent with our result for 
the zigzag edge. We speculate that the sample used in 
their experiment is similar to a nanoribbon with an arm- 
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chair edge. This speculation is reasonable because arm- 
chair edges are more frequently observed in experiments 
than zigzag edges. 28 ) Casiraghi et al. performed Raman 
spectroscopy on graphene edges and observed a small 
rcdshift of the G peak near the edge accompanied by a 
decrease in the linewidth of the G peak. 29 ) This behavior 
of the G peak is consistent with that of the zigzag edge 
since it is only the TO mode without broadening (which 
is related to the imaginary part of the self-energy) that 
can be Raman active. 

The cutoff energy E c appearing in eq. (26) may be 
determined from a tight-binding lattice model. For peri- 
odic graphene, by taking into account the contribution 
of all the possible electron-hole intermediate states in the 
Brillouin zone, we can have E c rts 770 (20 eV), which is 
larger than the value adopted in eq. (26). The value of 
E c for a graphene sample with an edge may be differ- 
ent from that for a periodic graphene sample without an 
edge. In fact, for a nanoribbon, a tight-binding calcula- 
tion 25 ) shows that the energy difference between the LO 
and TO modes is approximately 30 cm" 1 , which corre- 
sponds to E c ~ 6 eV. Thus the value of E c depends on 
the geometry of the system. Since we have considered a 
large graphene sample with an edge, we assumed that an 
appropriate value of E c is between 6 and 20 eV, and we 
chose 10 eV, which is of the order of half of the tt band- 
width. Because E c is not an experimentally controllable 
parameter, we consider that, in order to verify our re- 
sults, it is essential to observe the Ep dependence of the 
G band spectra near the edge. 

In conclusion, the el-ph matrix elements for the Raman 
intensity and Kohn anomaly near the edge of graphene 
were derived by adiabatic calculation, and then pertur- 
bation treatment was applied to the nonadiabatic parts 
of the phonon self-energies. The zigzag edge causes in- 
travalley scattering and the y-componcnt of the pseu- 
dospin vanishes (a y ) = for the standing wave. The 
Raman intensity of the LO mode and the Kohn anomaly 
of the TO mode are negligible owing to (a y ) = 0. On the 
other hand, the armchair edge causes inter valley scat- 
tering and the pseudospin does not change its direction. 
However, owing to the interference between two valleys 
originating from the el-ph interaction, the Raman inten- 
sity and Kohn anomaly are negligible only for the TO 
mode. The Raman intensity is enhanced when the po- 
larization of the incident laser is parallel (perpendicular) 
to the armchair (zigzag) edge. The difference in the be- 
havior of the pseudospin with respect to the zigzag and 
armchair edges is the origin of the asymmetry between 
the LO and TO modes. Our results are summarized in 
Table I. 
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Table I. Dependences of the Raman intensities and Kohn anoma- 
lies on the r point optical phonon modes. The symbols O an d 
X for Raman intensity and the Kohn anomaly represent 'oc- 
currence' and 'absence', respectively. There is asymmetry be- 
tween the Raman intensity and Kohn anomaly, that is, the Kohn 
anomaly occurs only for the LO mode, while the mode with a 
strong Raman intensity changes according to the edge shape. 
Raman intensity is enhanced when the polarization of the inci- 
dent laser light is parallel (LO) to the armchair edge or when it 
is perpendicular (TO) to the zigzag edge. 



Position 


Mode 


Raman 


Kohn 


Polarization 


zigzag 


LO 


X 


O 


X 




TO 


O 


X 


O 


armchair 


LO 


O 


O 


O 




TO 


X 


X 


X 


bulk 


LO 













TO 












Appendix A 



of Edge States to 



Correction 
eq. (10) 

we exactly calculate the x-componcnt of the 



Here 

pseudospin, (^kl^l^k)' with k y 
wc have 



0. When k v 



0. 



cos0(k). (A-l) 

It should be noted that this expression holds for extended 
states. The states with k v = are divided into two states, 
extended states and edge states, 21 -* depending on the sign 
of fca;. 30 ) For the K point, the edge states satisfy k x < 0, 
while the extended states satisfy k x > 0. Since the edge 
states are pseudospin polarization states, that is, they 
are eigenstates of a z , then the matrix element of a x with 
respect to the edge states vanishes. Thus, the exact form 
is given by 

cos6»(k) (fcyT^O), 
sign(fc x )cos#(k) (k y = 0), 

where sign(fc x ) = 1 when k x > and zero otherwise. By 
neglecting this complication, we obtain eq. (10). 



(A-2) 



Appendix B: Derivation of eq. (25) 

In this section, we derive eq. (25). 
First, using eqs. (17) and (22), we obtain 

(*kl^Lo m |*k> - -«^sm0(k). (B-l) 

Next, we consider the real part of the self-energy by set- 
ting S = in eq. (23). At zero temperature, we can set 
/h - / e = 1 for > |-E F |, otherwise / h - f e = 0. Then, 
the self-energy of the LO mode is written as 



Be[n(cj,E F )}=2(v F A«) 2 Y^ 



sin' 



0(k) 



sin' 



0(k) 



(B-2) 

where Y^' indicates that the summation is taken over 
states satisfying > \Ep\. Since the y-axis (a;- axis) is 
parallel (perpendicular) to the armchair edge, we use a 
periodic boundary condition for k y and an open bound- 
ary condition for k x . Then we have k y = 2nn y /L y and 
k x — nn x /L x . The summation over possible electron- hole 



8 J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



pairs can be rewritten as 



— - I dk x / dky 

2tH J Q J_ kc 



(B-3) 



where V = L x L y , k c is the cutoff momentum, and the 
factor of 2 originates from the degeneracy with respect to 
the K and K' points. Substituting eq. (B-3) into cq. (B-2), 
we have 

Re [n(w, E F )} = -J- J sin 2 6d6 



EdE { — 



1 



Hlo + 2E 



(B-4) 

where we have changed the integration variables from 
(k x ,k y ) to {E,k) by using k x ~ kcos9, k y = fcsin#, 



and E = hvpk. Using J^^ 2 sin 2 @d6 = ir/2 and J x/(x + 
a)dx = x — ahi\x + a\, we obtain eq. (25) when E c Hui. 

To calculate A%, we have used = g gu(uj)/a cc 
where a cc = 1.42 A. Here g g is the off-site el-ph ma- 
trix element and u(u>) is the amplitude of the phonon 
mode. We adopt g s = 6.4 eV. 31 ) A similar value is 
obtained by a first-principles calculation with the local 
density approximation. 32 ' We use a harmonic oscillator 
model which gives u(lj) = ^h/2M c N u uj, where M c is the 
mass of a carbon atom. Using hw = 0.2 eV, we obtain 

y/TQJQ/ti\ « 2 x lO" 2 !" 1 . 
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